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The structural control of silicon nanocrystals is an important technological problem. Typically a 
distribution of nanocrystal sizes and shapes emerges under the uncontrolled aggregation of smaller 
clusters. The aim of this computational study is to investigate the evolution of the nanocrystal 
electronic states and their optical properties throughout their aggregation stages. To realistically 
tackle such systems, an atomistic electronic structure tool is required that can accommodate about 
tens of thousand nanocrystal and embedding lattice atoms with very irregular shapes. For this 
purpose, a computationally-efficient pseudopotential-based electronic structure tool is developed 
that can handle realistic nanostructures based on the expansion of the wavefunction of the aggregate 
in terms of bulk Bloch bands of the constituent semiconductors. With this tool, the evolution of 
the electronic states as well as the polarization-dependent absorption spectra correlated with the 
oscillator strengths over their aggregation stages are traced. The low-lying aggregate nanocrystal 
states develop binding and anti-binding counterparts of the isolated states. Such information may 
become instrumental with the maturity of the controlled aggregation of these nanocrystals. 
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I. INTRODUCTION 

Embedded semiconductor nanocrystals (NCs) can be 
fabricated by several different techniques such as ion im- 
plantation [l|, [4I , sputtering [3|] , plasma enhanced chem- 
ical vapor deposition [3]. However, their structural con- 
trol is still an important technological problem. Dur- 
ing the annealing process, the excess silicon atoms start 
forming clusters and with their gradual aggregation into 
nanocrystals, a distribution of sizes and shapes emerges. 
The annealing temperature and duration determine the 
degree of the cluster aggregation and eventually the av- 
erage NC size distribution. Being directly controlled by 
this distribution 4, 5] the luminescence properties are 
intimately affected by this process of NC aggregation 
driven by atomistic nucleation and growth kinetics [6|. 
It has been generally demonstrated that the increase in 
the Si concentration of the substoichiometric silicon ox- 
ide films, as well as the annealing temperature lead to 
larger NCs and hence to the red shift of their lumines- 
cence spectra [J]. However, the experimental studies are 
hampered by the large size distribution existing in typ- 
ical samples. The soul of the present work is based on 
the controlled agglomeration of Si NCs. This is currently 
not achievable by the PECVD, sputtering or ion implan- 
tation techniques all of which are kinetically controlled. 
One of the most sophisticated fabrication tools for this 
purpose is the extreme ultraviolet interference lithogra- 
phy [7] where silicon-based two-dimensional quantum dot 
arrays with nanometer size control have been successfully 
demonstrated. 

In this work we tackle the electronic structure and 
the optical properties of Si NCs along their aggregation 
stages from a computational point of view. For a realis- 



tic account of the aggregation process, an atomistic elec- 
tronic structure tool is required that can accommodate 
about tens of thousand nanocrystal and embedding lat- 
tice atoms with very irregular shapes. However, currently 
highly sophisticated ab initio density functional theory- 
based plane wave approaches can only tackle systems be- 
low thousand atoms even with the help of ultrasoft pseu- 
dopotentials 8]. For this purpose, a computationally- 
efficient pseudopotential-based electronic structure tool 
is developed that can handle realistic nanostructures 
based on the expansion of the wavefunction of the ag- 
gregate in terms of bulk Bloch bands of the constituent 
semiconductors. With this tool, the evolution of the elec- 
tronic states as well as the absorption spectra and in par- 
ticular the oscillator strengths of transitions of interest 
over their aggregation stages are traced. 



II. THEORETICAL DETAILS 

We employ a computationally-efficient 

pseudopotential-based electronic structure tool that 
can tackle on the order of tens of thousand atom 
nanostructures based on the so-called linear combination 
bulk bands recipe of Wang and Zunger [^, [lOl- The 
crux of the technique is to expand the wavefunction of 
the NC in terms of bulk Bloch bands of the constituent 
semiconductors ; 
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where, vF -(f) is the cell-periodic part of the Bloch states, 
N is the number of primitive cells within the computa- 
tional supercell, C"' - is the expansion coefficient set to 

n,k 

be determined and a is the constituent bulk material la- 
bel which points to either the NC core or the embedding 



media. The atomistic NC Hamiltonian is given by 

- + 

a.Rj.a 



inary part of linear dielectric function becomes 
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where m is the free-electron mass, W^{Rj) is the weight 
function that usually takes values between to 1 depend- 
ing on the type of atom at the position Rj + d^ (which 
becomes instrumental especially at the interfaces) and v"^ 
is the screened spherical pseudopotential of atom a of the 
material a. We use semiempirical pseu dop otentials for Si 
and Ge taken respectively from Refs. |ll| and [1^. This 
is one other advantage of this simpler approach over the 
more accurate ah initio techniques which run into well- 
known band gap shrinkage due to local density approxi- 
mation ^] restoration of which is not straightforward. 
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FIG. 1: Empirical pseudopotential band structure of the Si 
(solid) and the wide band gap matrix (circles). 

An important issue is the choice of the host matrix 
material. If the NC is surrounded by vacuum, this corre- 
sponds to the free-standing case. However, the dangling 
bonds of the surface Si atoms lead to quite a large 
number of interface states which adversely contaminate 
the effective band gap region of the NC. In reality Si 
NCs are always embedded into a wide band gap host 
matrix which is usually silica [3, 0]- For this reason, we 
embed the Si NCs into a wide band gap medium having 
proper band alignment so that both electrons and holes 
of Si will experience quantum confinement. This is 
illustrated in Fig. [1] showing the bulk band structure 
of Si and the wide band gap medium with a band gap 
above 6 eV. 

With the possession of the electronic wavefunctions of 
the NCs, their linear optical properties can be readily 
computed. Within the independent-particle approxima- 
tion and the artificial supercell framework [13j the imag- 
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where, i = x,y,z denotes the cartesian components of 
the dielectric tensor and 
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is the oscillator strength of the transition. In these ex- 
pressions V is the volume of the supercell and Nk is 
the number of k points chosen within the Brillouin zone, 
the label v (c) corresponds to occupied (empty) valence 
(conduction) band states. Note that Im{eii} implicitly 
depends on the filling factor / — Vnc/^ which is the 
ratio of the NC and the supercell volumes. To obtain 
a more general quantity an effective dielectric function 
can be introduced as Im{eii}// [lj|. This is the form 
we prefer to present our results in the next section; for 
any specific filling factor one can readily infer the corre- 
sponding dielectric function. As another technical detail 
a Lorcntzian broadening of 0.2 cV of the Dirac delta func- 
tion in Eq. ^ is used. Finally, the absorption coefficient 
a{uj) is related to the imaginary part of the dielectric 
function through |15| 



Ime(a;) 



TlrCtQ 



LO 



a{huj) 



(5) 



This expression is in SI units (unlike the previous expres- 
sions) where Ur is the index of refraction, c is the speed 
of light and eo is the permittivity of free space. 

Even though there is no geometrical restriction of the 
approach, we consider the aggregation of two equal- 
diameter spherical NCs as illustrated in Fig. [21 Here, 
the isolated diameters are taken as Dg and the center-to- 
center distance between the two NCs as XgOg- Assum- 
ing no strain generation during the aggregation process 
(0 < a^s < 1), the diameter of each NC should expand 
from Dg to Da to preserve the overall volume and the 
number of atoms. A simple geometrical consideration in- 
dicates that the expansion factor y — Da/Dg should be 
determined by solving the following cubic equation 
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for a given Xs value that describes the degree of aggrega- 
tion. This results in some underestimation of the number 
of atoms due to those lost in the reduced surface area of 
the aggregate; therefore we slightly enlarge the aggre- 
gate size until we accommodate the correct total number 
of atoms. 
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FIG. 4: Evolution of the low-lying valence NC states; the Xs 
values are indicated above the spectra. 
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FIG. 2: Illustration of the aggregation stages of two equal- 
sized NCs 
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FIG. 3; Evolution of the low-lying conduction NC states; the 
Xs values are indicated above the spectra. The two states 
labelled as -a and -b are considered in Fig. [5] 



III. RESULTS 

Two NCs each with isolated diameters of 1.6 nm are 
assumed to approach together along a common x— axis. 
The evolution of conduction and valence NC states are 
shown in Figs. [3] and [3] respectively, with the aggregation 
being quantified by the x^ parameter. In these scales 
some of the degeneracies in both the conduction and va- 
lence NC states may not be resolved when spherical sym- 
metry is restored for Xg = and Xg — » oo. However, in 
intermediate aggregation stages these degeneracies are 
further lifted. Also note that the conduction NC states 
are mainly derived from 0.85X minima of the bulk Si, 
hence tend to be six-fold degenerate for the low-lying 
states. Mixing with other bulk band states again lifts 
the exact degeneracy. From both of these figures, the 
splitting of the isolated NC states into binding and anti- 
binding sets can be observed. This is illustrated in Fig.[S] 
where the wave functions of the two states marked as -a 
(for anti-binding) and -b (for binding) are plotted clearly 
justifying their identification. 

The effect of NC aggregation on the optical properties 
is shown in Figs. [H] and [7] choosing Xg — 0.8 and cases 
for comparison. We should remark that the embedding 
medium has no contribution within the energy window 
(below 3.6 eV) considered here and all absorption is due 
to real transitions between NC states over the effective 
band gap. The oscillator strengths are plotted for the xx 
component of the tensor where x is the direction along 
the aggregation. Since the aggregation of the two NCs 
leads to a larger NC, transitions are shifted to lower 
energies, as expected (cf. Fig.[S]). Figure [7] demonstrates 
that XX component of the imaginary part of the dielectric 
function is enhanced whereas the other two components 
attain a lower but equal value which is also the case for 
the larger and the fully aggregated and spherical Xg = 
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FIG. 5: Wave function amplitude square plots of the binding 
and anti-binding states indicated in Fig. |3]for the a;^ = 1 case. 
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FIG. 6: The xx component of the oscillator strength tensor 
of the low-lying transitions for Xs — 0.8 and cases. 



case. Note that the energy for the local peak of the xx 
component of the Xs — 0.8 absorption spectra correlates 
with strong oscillator strengths in Fig. [6] around 3.15 eV 
which are absent for Xs = 0. Obviously, in the case 
of random growth processes the macroscopic average 
will wash out such polarization effects, however for 
more controlled fabrication techniques may reveal these 
features. 



FIG. 7: The imaginary part of the effective dielectric tensor 
associated with the previous figure. 



IV. CONCLUSIONS 

A pseudopotential-based atomistic electronic structure 
tool is developed to realistically investigate the aggrega- 
tion process of NCs. With this tool, the evolution of the 
electronic states as well as the absorption spectra and in 
particular the oscillator strengths of transitions of inter- 
est over their aggregation stages are traced. The predom- 
inantly X-derived six-pack groupings of conduction NC 
states are observed. The low-lying aggregate NC states 
are identified as the binding and anti-binding splittings 
of the isolated states. The effect of the elongation of 
the NCs in the initial phases of a controlled aggregation 
is reflected as a polarization-dependent absorption spec- 
tra. Such information may become instrumental with 
the maturity of the controlled fabrication of these NC 
aggregates. 
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